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Abstract 

A minimal model for coherent transport through a donor/acceptor molecular junc- 
tion is presented. The two donor and acceptor sites are described by single levels 
energetically separated by an intramolecular tunnel barrier. In the limit of strong 
coupling to the electrodes a current rectification for different bias voltage polarities 
occurs. Contacts with recent experiments of molecular rectification are also given. 
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1 Introduction 



In recent years, molecular electronics has gained considerable attention due 
to the viability to measure quantum transport observables through single 
molecules. In fact, various effects typical of the physics of mesoscopic electron 
systems ranging from Coulomb blockade effects with weak coupling between 
the molecule and the leads to coherent transport through strongly bonded 
molecules have been observed [1,2]. However, it is desirable to obtain clear 
relationships between truly intrinsic molecular properties and measured ef- 
fects. For example, vibron assisted tunneling [3] or current switching due to 
a controUably induced conformational change of the molecular shape [4,5] 
have been observed. Another idea, firstly proposed theoretically by Aviram 
and Ratner [6], is to create molecules able to rectify electric currents, i.e. 
I{U) 7^ —I{—U). This can be either achieved by distinct asymmetric molecules 
(e.g. asymmetric tunneling barriers between the contacts and a central part 
of the molecule [7]) or by an intrinsic modification of the electronic properties 
of the molecules [8,9]. 
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2 System and Method 



Having in mind the experiments by Elbing et al. [9], where a molecule con- 
sisting of two broken vr-conjugated donor /acceptor subunits connected to gold 
electrodes via thiol groups showed a diode-like behavior, a minimal model for 
quantum transport through such a molecular system is presented. Though 
lacking ah initio accuracy, our model gives a basic insight to the mechanism 
of molecular rectification. 

We describe the donor-acceptor molecule by two single levels separated by a 
tunneling junction. The on-site energies are e^i and (in the case of the 
isolated molecule), respectively (see Fig. 1). These energies are measured rel- 
atively to the Fermi energies of the electrodes, which are both set arbitrarily 
to zero. The two parts of the molecule are coupled by a weakly conducting 
bridge, which can be formed by a a-bond [6] or an almost broken biphenylic 
vr-bond [9]. Under the external bias voltage [/, which is applied symmetrically 
on the two contacts, the chemical potentials of left and right lead are moved: 
/xl = — /iR = eU/2. Furthermore, the energy levels of the molecular sites are 
shifted as 

en {U) = e?i + r^^eU, £22 {U) = + T]<,dJ. (1) 

To implement the fact that the bottleneck tunneling barrier in the system 
would tend to pin these levels to the nearby electrode, we assume that one 
fourth of the applied voltage drops at the contacts between the molecule and 
the leads and one half of the bias voltage drops in the middle of the system 
between the two parts of the molecule at the intra-molecular barrier, i.e. ?7i = 
— ?72 = 0.25, this assumption is then lifted at the end of this paper when 
discussing different strengths of coupling between electrodes and molecule. 




(a) ?7 < (b) [/ = (c) [/ > 



Fig. 1. Sketch of the energy levels at negative (a), zero (b) and positive (c) bias 
voltage. The dotted hnes indicate the position of the Fermi energy. Since the molec- 
ular levels themselves are affected by the bias voltage, Eq. (1), the distance and 
positions of the levels are strongly sensitive to the bias polarity. 
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The Hamiltonian describing the molecule is given by 

a,/3=l,2 



(2) 



whereas the diagonal elements e^a are defined by the on-site energies Saa (U) 
and the off diagonal elements ei2 and 621 are given by the hopping parameter 
between the two levels, —7, which is a measure for the intcr-sitc coupling. In 
the two dimensional base of the localized orbital operators the Hamiltonian 
matrix simply writes: 



H 



M 



^ -7 £22 (U) J 



(3) 



The total Hamiltonian contains two additional terms due to the electrodes and 
their coupling to the molecule: H — Hu + -f^M-ieads + -f^ieads- To calculate the 
spectral and transport properties of the molecule coupled to leads, we use the 
nonequilibrium Green functions technique for finite bias voltage without in- 
teraction. The retarded Green function matrix of the molecular region dressed 
by the electrode self-energies reads [10] 
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are the self-energy matrices of left and right lead, which already lift the molec- 
ular resonances away from the real energy axis, surrogating the need of a 
small imaginary shift applied in the definition of the retarded Green function. 
For the self-energies we use the wide band approximation, which assumes a 
purely imaginary energy-independent self-energy. Al,r describe the hopping 
between the contacts and the two energy levels. Both Al and Ar are assumed 
to be constant positive real numbers with no energy dependence (wide band 
approximation). This is a reasonable assumption when thinking of gold elec- 
trodes whose bands are more extended than the molecular active energetic 
window. The transmission probability is obtained from the Green function by 
the Fisher-Lee relation T{E,U) = [Ti^G'T^G'^} [11,12] with the matrices 

Fl = i (^Sl ~ ^l) Fr = i ^Sj^ — Sf^^ as the anti-Hermitian parts of the 
self-energy matrices of the contacts. 
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In the case of the two site-model studied here, it is easy to obtain an analytical 
expression for the transmission probability as a function of voltage and charge 
injection energy: 



T{E,U) 



B 



A 



4AlAr =4AlAr7V(A + S) 

\e - en ([/)) {E - 622 (U)) - 7' - AlAr 
[Al {E - en (U)) + Ar - €22 {U))f . 



(6) 



Knowing the transmission probability, the current through the system can be 
obtained by the relation [12]: 



Here, /l r (E) = f {E — /xl,r) ai'e the Fermi functions of left and right elec- 
trode depending on the bias voltage and the temperature, which in our cal- 
culations is set to 30 K. To understand the evolution of the molecular levels 
under the influence of the applied bias voltage and the coupling to the leads, 
we analyze the density of states (DOS), especially the localized density of 
states (LDOS) projected on the two sites of the molecule: 



Additionally, we define Ados = LDOS1-LDOS2 and Hdos = LDOSi-LDOSa- 
Looking at Ados, one can find the localization of a molecular level: positive 
Ados implies that the state is more localized on site 1, negative Ados points 
to a pronounced state on site 2. IIdoSi on the contrary, enhances with the 
equipartition of molecular states on the two sites. 



3 Results and Discussion 

The case of two identical sites eJi = €22, independently of their value, leads 
straightforwardly to a symmetric behavior of the current with respect to an 
inversion of the applied bias, i.e. / (U) — —I {—U). First, we want to study the 
case of two levels with different energies coupled symmetrically to the leads 
(Al = Ar = Aq). In the upper panel of Fig. 2, the difference of the DOS 
projected on the two levels is shown as a function of the applied bias voltage 
and the energy: the two lines show the variation of the two molecular levels 



I{U) = ^jdET {E, U) (/l {E) - /r {E)) 



(7) 



LDOSi {E, U) = -^^ImTriG^ = 
27r 

LDOS2 {E, U) = -;^ImTr2G"^ = 




(8) 
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Fig. 2. Difference (upper panel) and product (lower panel) of the localized density 
of states of site 1 and site 2 as a function of the bias voltage and the energy. 
The white lines indicate the bias window. In order to emphasize the molecular 
states, a logarithmic nonlinear scale has been applied. (Parameters: e^i = —0.3 eV, 
e^2 = 0.6 eV, Aq = 0.02 eV, 7 = 0.3 eV, rji = -r]2 = 0.25). 

which are broadened by the couphng to the leads. The locahzation of the levels 
is indicated as follows: red resp. black color of the central part indicates that 
the level is situated on site 1, blue points to a level localized on site 2. For large 
negative bias voltages, the two levels are far apart. With increasing voltage, 
the two levels approach each other, but do not cross and move away from each 
other again. This repulsion is quantified by the inter-site coupling 27. At this 
voltage. Ados is zero whereas IIdos takes its maximum (lower panel). After 
the avoided crossing, a change in the localization of the energy levels takes 
place. 
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Fig. 3. Transmission function T {E, U) (upper panel, in a logarithmic scale) and 
current / (U) (lower panel) for a two-site system. The white lines in the upper 
panel indicate the bias window. Different rectification degrees are obtained for 
different asymmetry parameters a, with Al = Aq and Ar = aAg. (Parameters: 
-0.3 eV, e^2 = 0-6 eV, Aq = 0.02 eV, 7 = 0.3 eV, rji = -772 = 0.25). 
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Fig. 3 shows the calculated transmission function and the current flowing 
through the system (thick blue line in the lower panel) calculated using Eq. (7). 
The bias window, indicated in the plot of the transmission function by the 
two white lines, is determined by the integration interval fixed because of the 
Fermi functions at low temperature. A step in the current-voltage curve ap- 
pears, when a level enters the bias window, as it happens once at negative and 
once at positive bias voltages. If Aq is small compared to the level spacing, 
i.e. for relatively sharp molecular levels, the positions of the steps is given 
by the crossing of the eigenenergies (this can be determined analytically by a 
diagonalization of the Hamiltonian in Eq. (3)) and the lines defining the bias 
window. In our case, as depicted by the solid grey vertical lines in Fig. 3, this 
takes place at U = —1.41 V in negative bias direction and aX U = 0.56 V and 
U = 1.01 V for positive bias. The height of the steps however depends cru- 
cially on the distance between the two levels: in the region of negative voltage, 
the two states are far apart and the height of the steps is therefore relatively 
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small, whereas for positive bias voltage, the levels are much closer to each 
other, which results in a more resonant state and a much higher current. This 
explains the observed rectification in our model. The highest current is reached 
where the two levels are closest to each other. The peak in I{U) appears at 
the bias voltage {€22 — /e (r?i — 772) (in our case U — 1.80 V, depicted by 
the dashed grey line in Fig. 3) which can be obtained by a minimization of 
the energy gap between the two eigenenergies of the bare molecule [13]. For 
higher voltages, the distance between the levels grows again and the current 
decreases. 

So far, we only considered symmetric couphng of the molecule to the leads. 
However, in experimental investigations of electronic transport through sin- 
gle molecules, the strength of the bonds between the molecule and the leads 
is mostly not well defined. Therefore, it is important to study the case of 
the molecule coupled asymmetrically to the electrodes, i.e. Al = olAq and 
Ar = orAq. The thin lines in Fig. 3 represent the values of the current 
for asymmetric coupling ol = 1 and different values of Cr = a, resulting in 
Al = Aq and Ar = oAq. The other model parameters remain unchanged, 
especially rji = —r]2 = 0.25. In this case, the general shape of the current 
curve does not change; the steps and the maximum in the current voltage 
curve do not change their positions, only the absolute values of the current 
are different due to the varying coupling to the leads. To understand the cor- 
rect voltage drop along the molecule, one should use the Poisson equation with 
boundary conditions given by the applied bias voltage. This goes beyond the 
scope of this paper which concentrates on a simplified picture of the molecular 
electronic structure. Thus, in order to improve the plain effect of asymmetric 
coupling to the shape of the voltage drop along the system, we assume a larger 
voltage drop at the weaker contact. This changes as a consequence the values 
of the factors 771^2 defined in Eq. (1) according to, e.g,: r/i = (3aR — o-l) /8 
and ri2 = (or — 3ai^) /8, which for the case = or = 1 would restore the 
previously used values of 771^2- The results for this refined model are shown 
in Fig 4, the thick blue line again representing the symmetric case. Here, the 
position of the steps is shifted to different voltages. As it becomes evident 
from a comparison of Figs. 3 and 4, different shapes of the voltage drop imply 
a shift of the position where the current significantly changes, i.e. where the 
steps appear. Rectification effects due to the asymmetric coupling as shown in 
Fig. 3, have been already reported in the literature [14] but we here we want 
to stress that the precise value of the current jumps are not a mere electronic 
structure effect, but also depend on the profile of the electric field under the 
fixed applied bias voltage. 
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Fig. 4. Current voltage characteristics for two cases of asymmetric coupling. In 
the upper panel, ol = 1, the lower panel shows the case or = 1. (Parameters: 
= -0.3 eV, e^^ = 0-6 eV, Aq = 0.02 eV, 7 = 0.3 eV, r?i = -772 = 0.25). 

4 Conclusions and Outlook 



In this paper, we have shown a minimal model for rectification effects in 
coherent transport through single molecules which was inspired by recently 
published experiments [9]. An extension to this model includes charging ef- 
fects [15], taking into account on-site correlations on the two molecular sites, 
which becomes important in the case of a weak coupling between the molecule 
and the electrodes. These implementations though do not essentially change 
the nature of the effects observed in the picture discussed here, except for a 
natural splitting of conductance steps due to the lifting of the spin degeneracy 
induced by the correlation effects. Additionally, density functional theory can 
be used to calculate the positions of the energy levels of the molecule and their 
exact behavior under each bias voltage applied to the molecule, combined with 
a calculation of the correct voltage drop distribution at this voltage solving the 
Poisson equation. Together with nonequilibrium Green functions, this proce- 
dure allows to calculate the transmission T {E, U) and the current I (U), and 
not only the electronic structure under an applied voltage, in a truly ab initio 
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way. Such study is currently under investigation [16]. 
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